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Abstract 

We tackle the problem of the accelerating universe by reconsidering the most general form of the metric 
when the speed of light is allowed to evolve with time in a homogeneous and isotropic universe. A new varying 
speed of light (VSL) model naturally emerges. We find it unconvenient to perform a general coordinate 
transformation to obtain the usual constant gtt as such an operation would mask the relation existing 
between the evolution of the scale factor and the frequency shifts. In the model proposed the expansion rate 
and the acceleration of the universe turn out to be apparent effects induced by the evolution of the speed of 
light. The model is beneficial in that no sort of exotic (and so far unobserved) fluids, not even a cosmological 
constant, are needed for our solutions to be compatible with observations. Only dust and radiation are put 
into the energy-momentum tensor and these are found to be sufficient to reach the critical density for the 
model and therefore to obtain a spatially-flat universe. The field equations for the model are derived and 
solved. Among others, one fashinating possibility is that of an eternally bouncing universe. The relation 
with the varying-Q EM results is discussed. 
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1 Introduction 



"It conflicts with one's scientific understanding to conceive a thing which acts but cannot be acted 
upon. " 

A. Einstein on the concept of absolute space. 

Many authors have recently proposed cosmological models where a variation with cosmological time of the 
universal constants, in particular of the speed of light, is hypothesized as a viable alternative to the widely 
accepted inflationary scheme in order to solve the classical problems of cosmology (see pQ and references therein 
for a comprehensive review of the topic) . 

Attempts in this direction started about ten years ago with Moffat [2j who considered a sudden change in 
the speed of light and discussed the cosmological implication of such a phase transition. The same kind of 
variation was used a few years later by Albrecht and Magueijo to propose their c-varying model that, under 
certain conditions, presents no cosmological problems. A model for the universe with c evolving according to 
a power law was proposed in the same year by Barrow [4], who also looked for solutions to the field equations 
and set the constraints under which his solutions solve the classical problems. 

Other authors also opted for a continuously varying c(t), see for instance [5], where Belinchon and Alfonso- 
Faus introduced the Planck's constant in the field equations via the Stefan-Boltzmann law, discussed the gauge 
invariance for the Schrodinger, Klein-Gordon-Fock and Maxwell's equations and showed that their model does 
not present the so-called Planck's problem. See also the "revised" version [6]. In other works (see e.g. Bel in 
[Zj), the speed of light directly acquires the role of the reciprocal of the scale factor, with the metric divided by 
c 2 , and the dimensionality of a length squared preserved through multiplication by an appropriate constant. 

Another interesting way to introduce a varying speed of light is the bimetric theory of gravity proposed 
by Clayton and Moffat [5] where the issue of diffeomorphism invariance is elegantly approached by assuming 
the presence of two different metrics which implies a different speed for gravitons and photons. This model is 
capable to solve the classical cosmological problems as easily as inflation does but, contrary to inflation, it does 
not need fine tunings or unnaturally flat potential^]. More recently [9], Clayton and Moffat have shown that 
their bimetric theory predicts a scale-invariant spectrum for cosmological perturbations which is compatible 
with CMB data but distinguishable from the inflationary prediction. 

Despite some evidences for a varying fine-structure constant a EM [12} [131 El US] (contested for instance 
in [TB] and rebutted in p2]), which at the time being represent the only observative evidence for the variation 
with redshift (and therefore with time) of (at least) one "universal constant", there are some authors who contest 
the very meaning of stating that the speed of light might be not constant (see e.g. in [18]). The most popular 
objection seems to be that, after all, it all boils down to the choice of a system of units through which one can 
see a varying-a theory as if h or e were time dependent instead of c, or together with c, or finally in any possible 
combination of the three (see also [19] for a reactionary view on varying-c theories). This objection deserves 
consideration even if, in the general case, the appropriate redefinition of units has to be continuous in, say, the 
cosmic time. 

After some 15 years from their "discovery", the discussions about the meaning of VSL theories is going on. 
A recent opposers' viewpoint may be found in [20] where a sort of general check list aimed to decide if a VSL 
model merits consideration at all is proposed. This list however does not seem to take into account the different 
ways and purposes to introduce a VSL mechanism, and the true intention seems to show that introducing a 
VSL is equivalent to upset the very basics of modern physics, which is certainly not the case in general and 
surely is not true for the model presented in this paper. 

The arguments presented in [20] have been quickly replied in [21] where a particularly illuminating example 
can be found (see section II). It relates to acceleration of gravity, the once-believed "constant" g, and makes clear 
how the arguments put forth by the opposers of VSL theories could have well applied at the time of Newton 

1 For a review about the most important problems with the inflationary paradigm see, for instance, |10j . 
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against any proposals that g were not constant. This wonderful example should make clear to anyone how such 
choices really depend on the knowledge of physics at a given time. Saying that c is constant because so we 
have decided (so we have defined units) is not different from insisting that g is constant because so Galileo had 
decided. As Galileo was wrong because at his time he could not measure g at a large distance from the Earth 
surface, so we may well be wrong because in no way we can know now what the speed of photons was 10 9 years 
ago (see the next section for what we exactly mean with "the speed of photons has varied") . We do not know if 
the speed of light has varied or not, but it is sure that it is a possibility that worths be investigated. 

As said, the reported evolution of a EM may also be attributed to a varaition of e or h. An example of a 
time-dependent a EM described by a varying-e theory, we address the works by Bekenstein |221 [23] and their 
generalization by Barrow, Magueijo and Sandvik [24] (the BSBM model). Sometimes the variation of a EM has 
been considered together with that of G. Time-dependent G theories date back to the well-known works by 
Dirac [25] and by Brans and Dicke [26]. More recently Barrow and Magueijo have studied the behaviour of 
Brans-Dicke cosmologies with a varying speed of light [27] . See also [28] for a parallelism between Brans-Dicke 
cosmologies and varying speed of light theories. 

While, for the reason given above, a variable fine-structure constant - whose variation has been recently 
contested, see [29] and references therein - cannot be taken as an evidence for a variation of the speed of light 
with cosmological time, it is certainly incorrect to consider a variable speed of light meaningless or tautological. 
The whys, clearly explained in [TJ (chapter two), are easily summarized as follows. Since a EM is a dimensionless 
quantity, no matter which system of units one chooses to adopt, if a EM evolves with time, it will do so whatever 
the units employed. As a consequence, if a EM = e 2 /ch depends on cosmic time then at least one of the universal 
"constants" c, h or e must depend on cosmic time as well. Suppose now to have a varying-a EM theory with, say, 
a time-dependent c, and that such a theory has a simple dynamical structure. It is then certainly true that, by 
means of an appropriate units transformation, one can always rewrite that varying-a EM theory so as to have a 
constant c and, say, an evolving e, but by doing so one may end up with an unnecessarily contrived formalism 
which was much more elegant and simple in the varying-c scheme. Besides, given that a EM is dimensionless 
and therefore cannot be made constant by means of any transformation of units whatsoever, it is clear that in 
the new system of units e and/or h would be time dependent, and it is not clear why the variation of e and h 
should be regarded as less problematic than the variation of c. 

I wish to remark that this "simplicity" argument is not (as some might think) too weak to justify a varying-c 
cosmology. Let me remind that even the most primitive definitions adopted in physics are justified by similar 
arguments, the most outstanding example being perhaps the difinition of time (see [30j par 1.5 page 23). We 
define time by means of regular and periodic phenomena not because it could not be done otherwise, but because 
we want a unit of time which is the same now and tomorrow, a feature which facilitates a comparison between 
physical phenomena which are distant in time and guarantees that the equations of dynamics look simple. 

It is also to be considered that the experimental tests proposed in [31] show how varying-c and varying-e 
models may behave differently and could then in principle be distinguished if the experimental precisions were 
to reach a sufficiently high level. Proposed tests span from the weak equivalence principle to the dark matter 
needed to fit quasar data to earth-based gravitational redshift experiments (Pound- Rebka-Snyder) . For all these 
tests, Magueijo, Barrow and Sandvik [31J have shown that varying-c and varying-e models predict different when 
not opposite results. Finally, when a more complete and fundamental theory of physical phenomena will be 
available, we could eventually discover that one or more of what we now believe to be fundamental constants 
may in fact be derived quantities or just one of the factors combining to form a more fundamental constant. 
This, of course, would spoil that quantity of the role of "fundamental constant". 

To the best of our knowledge, assuming that the evidence for a variation of a EM will receive independent 
confirmations, it is presently impossible to predict which of the three constant a EM is made of is fundamental 
or not. We can only make guesses based on simplicity arguments. This principle of simplicity is the one that 
we will follow in the present work. It will be seen that one cannot rewrite this model in terms of a constant c(t) 
without paying the price of complicating the dynamics (e.g. by working with a continuously changing system 
of units) and making comparison with observations unnecessarily indirect and complicated. 

So far, works on VSL theories have mainly focused on showing that these theories are a viable alternative 
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to the inflationary scheme, especially when dealing with the solution to the classical problems of cosmology. In 
this paper we present a VSL model that provides an interesting solution to the puzzle of the acceleration of the 
universal expansion discovered some years ago p2j (for a more recent paper see [33j). 

Instead of introducing exotic and so-far-unobserved forms of dark matter-energy, we introduce only radiation 
or dust as sources in the energy-momentum tensor. We show that a universe filled with such "normal" content 
(and where the speed of light is varying) appears to accelerate even if the actual state of motion is not that of an 
accelerated expansion. Naively, this is due to the fact that instruments measure frequencies but the dynamics 
of the scale factor is "represented" by the evolution of wavelengths (see appendix A). Because of the relation 
\v = c, if the speed of light varies with cosmological time, the inverse proportionality between wavelengths and 
frequencies is lost and apparent effects arise. 

The paper is organized as follows. In section [2] we derive the general metric to be used under the approach 
proposed in the paper. In section[3]we introduce the concept of apparent scale factor and discuss how it is linked 
to the classical scale factor. In section 0] we give the expressions for the main tensorial quantities deriving from 
the metric proposed. Section [5] contains the phenomenological core of the paper: we show that the acceleration 
of the universe may be an apparent effect and derive under which circumstances such an effect takes place. In 
section [6] we solve the field equations in the case of spatial flatness and show that a universe composed of dust 
and/or radiation only actually meets the condition derived in section [5] under which an apparent acceleration is 
observed. In section [7] we perform some consistency checks to verify that the model proposed is both internally 
consistent and consistent with the relevant observative results. The summary is given in section [8j 



2 Postulates of the model and natural emergence of a VSL 

Most of the non-bimetric VSL models have been developed postulating that the speed of light has varied with 
cosmological time, either in a sharp phase transition or continuously. Such variation is assumed to affect the 
local conservation equation (the effect is "special relativistic") , but not the field equations [3J. The dependence 
of the speed of light c on cosmological time is simply accounted for by making the direct substitution c — > c(t) 
wherever appropriate. 

On the contrary, in our approach the form of the metric is affected by the evolution of c(t). As we have 
found that this position is considered questionable by many who are accostumed with standard cosmology, 
we dedicate this section to discuss the issue. As is well known, general relativity exhibits an invariance under 
general coordinate transformations. This implies, for instance, that we could in principle operate an appropriate 
transformation of the time variable so as to obtain a metric where c is constant. The advantage of operating 
such a transformation would clearly be the reduction of the degrees of freedom from two to just one, the scale 
factor a(t), and the usual FLRW metric would then be recovered. However obvious the naturalness of such 
an operation could appear, performing it is not at all compulsory, again because of general covariance. The 
choice is one of convenience and is dictated by the advantages and disadvantages implied by it. The question 
should be: does such an operation imply a global advantage or a global disadvantage? The answer depends of 
course on what the purposes of one's work are. We will find that, as far as comparison between the results of 
the model with observative evidences is concerned, operating the mentioned transformation would create more 
interpretative problems than the computational problems it would save us from. Therefore we find it simply not 
convenient to change the coordinate set, just as we do not use a meter stick whose length oscillates, increasing, 
for instance, during the day and decreasing during the night: while perfectly possible, comparison of lengths of 
other objects and its use in general would become an awkward task. 

Let us see examine the details of the problem. We start as usual with a four dimensional manifold endowed 
with some arbitrary metric and label a generic set of coordinates with x° , x 1 , x 2 and x 3 , all having the dimension 
of a length. The last three are chosen to be some "spatial" coordinates of convenience, for instance the spherical 
coordinates, while the a; coordinate must represent time in some way. In view of the cosmological application 
we next consider the restrictions imposed on it by assuming the validity of the cosmological principle. In the 
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case of spatial flatnes^l, the most general line element compatible with the cosmological principle is: 

ds 2 = b 2 (x a ){dx ) 2 - a(x°) 2 dr 2 (1) 

where b(x°) and a(x°) are some function of x° . Note that we have chosen the x° direction in such a way to make 
the gu elements vanis Here we may set b(x ) = 1 with no loss of generality as such an operation amounts to 
a redefinition of the meter, which cannot affect physical results (what we can measure is the ratio a(x°)/b(x )). 
We then obtain: 

ds 2 = (dx ) 2 - a{x°) 2 dr (2) 

where we have not bothered to rename the variable x° and the function a(x°). The choice made implies that no 
further simplifications can arise from changing units after this stage, otherwise one should put b(x°) or some 
other appropriate function back into the metric. 

The next step is to introduce time in the metric, both for sake of visualization (after all we do think in 
terms of time) and to make a comparison with the observation, in particular, with the redshifts or frequency 
shifts. What is the most convenient choice in our case, that is, in the case where one wants to consider the 
possibility of a time-dependent speed of light? One constraint is given by the fact that experience shows that 
physics is special-relativistic at a local level so that x a = tc must hold locally. Here c is the localljEI measured 
speed of light and t is the proper time along fluid worldlines (this time variable is global if space-time is flat) . 
The general transformation which is compatible with both this local limit and the cosmological principle is 
therefore x°(t) = tc(t) where c(t) reduces to the locally measured speed of light c. We remark that here we are 
substituting Lorentz invariance with local Lorentz invariance, see [33] . 

From x°(t) = tc(t) it follows that 

dx° = (c(i) + tc(t))dt (3) 
so that, in terms of the familiar coordinates, gu becomes: 

g tt = (c + tc) 2 (4) 

where we have dropped the time dependence stipulating that, from now on, c is always to be considered 
time-dependent unless otherwise stated. This implies that, under the restrictions imposed by the cosmological 
principle, the most general line element to be used in this varying-c model is 

ds 2 = (c + tcfdt 2 - a 2 dr (5) 

where, we repeat, spatial flatness has been assumed. 

At this point one would be tempted to redefine the time variable so as to obtain a constant speed of light: 
c(t) = Co. However, contrary to what happens in standard cosmology, performing such an operation happens 
to be not convenient in our case, and we dedicate the remainder of the paragraph to see why. 

We first discuss the importance of considering the local limit. By construction, the t variable used in the 
previous expression is the time of special relativity and, therefore, the time ticked by clocks at rest in a given 
reference frame. It follows that, in a cosmological setting, this time variable can be used as "the proper time in 
the comoving frame" that is, the time measured by clocks at rest with respect to the reference frame comoving 
with the cosmological fluid. If we assume that the relevant emission and absorption properties have not changed 
with time, then the use of such time variable guarantees that emission and absorption frequencies are numerically 
the same now and at any time. This property is important because it is precisely the requirement a time variable 
must fulfill for the observed redshifts to be completely ascribed to the cosmological expansion. 

We remark that this property remains true if and only if we use, up to a constant factor, the time variable 
given above. Indeed any coordinate change t — ► t' aimed to reduce gtn' to a constant would have to be time- 
dependent, t' = t'{t) (any space dependence is forbidden by the cosmological principle), and this would imply, 

2 The extension to the spatially curved case is straightforward. 

3 This is not different from what one does in standard cosmology upon choosing synchronous time as a time coordinate of 
convenience. 

4 With "locally" we are meaning here at a given x°, that is, on a given hypersurface of simultaneity. 
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among other things, that all emission frequencies, if expressed in terms of the new t' variable, could be different 
at different times t' , i.e. moving from one spacelike hypersurface to another. On the contrary, by using the time t, 
we are sure that we are using physical frequencies. It is evident that the crucial drawback of some transformation 
of the time variable would be that we could not anymore compare properties of radiation emitted in the past 
with what is known from earth-laboratory physics in a straightforward way. In particular any estimation of 
redshifts would necessitate a distance-dependent correction before comparison to earth-laboratory physics could 
be made. 

To see this point more clearly, let us attempt the conversion to an appropriate coordinate set to obtain a 
constant speed of light. We could for instance transform the time variable t in the following way: 

(c + ct old fdt 2 old — > cgd& w . (6) 

In this new reference frame the speed of light is constant but what are the consequences of such a transfor- 
mation? Probably the most evident is that while any given physical phenomenon (say, an atomic transition) 
which in terms of the i Q id variable lasted t \ d seconds at some epoch of the past will, by construction, also last 
fold seconds today, the same cannot be true in the i new frame or in terms of any other time variable different 
form t \d (up to a constant factor). The unpleasant consequences of the choice of some given i new are evident: 
the measured frequency shift experienced by any given radiation would have to be ascribed to two different 
contributions, one being the usual cosmological expansion effect and the other being the continuously changing 
emission/absorption frequency of the related atomic transition (due, as some put it, to the "evolving time unit") 
to which we compare the radiation to get the redshift. Only the first effect could be related to the dynamics 
of the scale factor and, as the measured redshift - which is measured by antennae which, of course, measure 
frequencies and not wavelengths - would carry the imprints of both effects, we would have to disentangle these 
two effects. Depending on the exact expression for i new = t new (t Q ] d ), the comparison with the observations would 
require a distance-dependent correction which could turn out to be an arbitrarily complicated and useless task 
as, because of general covariance, the results concerning the physical observables would in the end be the same. 

One conclude then that while in the i id frame the comparison between the theory and the observations 
is, by construction, straightforward, in the t new frame it could turn out to be unpredictably complicated. 
Since, because of general covariance, physical results, and in particular the evolution obtained for the physical 
observables, would turn out to be the same, there is no good reason to introduce a new, unconvenient time 
variable. This is why we will stick to the most simple choice and will keep working with the t id variable. 

This discussion should have made clear the simplicity argument we have been referring to: we use the 
varying speed of light frame and the related time unit because it is in this frame and only in this frame that 
the observed wavelength shift is purely cosmological. We close by noting that this remark allows to define what 
we do mean with a "varying speed of light" in a somewhat explicit form. Consider a physical phenomenon by 
which radiation of a given frequency is emitted. Think of the same phenomenon taking place at some time in 
the past and today. We have adopted the unit/frame by which the frequencies in the two cases are the same. 
Observe, for the two cases separately, the length a light ray can travel in, say, a unit of the time variable we 
have adopted. We will say that the speed of light has varied if that length is different for the two situtations. 



3 The apparent scale factor 

The most enlightening way to study the consequences of the modification to the tt coefficient of the metric 
induced by the variable speed of light and obtained in the last section is probably that of finding a relation to 
link the evolution of the scale factor a(t) (which represents the actual dynamics of the universe) with another 
variable, which conveniently represents the evolution of the scale factor as seen by an observer who is not taking 
the variation of the speed of light into account (henceforth referred to as the "constant-c observer") . We will 
denote such variable with A(t) and will refer to it as the apparent scale factor. 

As a first step toward this goal we have to re-express the evolution of a(x°) in terms of that of a(t). In fact 
the former will be seen to play the role of a bridge which connects the actual and the apparent scale factors. 
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From the evolution of a(x°) we will then be able to derive the evolution of the apparent scale factor. 

As already said, A(t) represents the evolution of the scale factor as seen by constant-c observers. In practice, 
this means that A(t) is what is usually obtained when estimating the evolution of a(t) from observations. As 
will be made clear later, it is the confusion between a(t) and A(t) which gives rise to the apparent effect of the 
acceleration in the expansion of the universe: A(t) "accelerates" while a(t) (the universe) does not. 

A preliminary remark will be helpful in understanding what follows. We recall that, whatever the method we 
may choose to perform any measurement on a cosmological scale, we have to use distant objects or, better, the 
radiation we receive from them. This radiation reaches us by travelling through space-time, along the worldline 
which separates the two events: (A) radiation emitted by the object, (B) radiation received by our instruments. 
Photons from these objects or from the LSS (Last Scattering Surface) have travelled their long way through 
space-time to reach our satellites and, in doing so they have experienced, event after event, the dynamics of the 
gravitational field, both the spatial and time metric coefficients continuously varying along their paths. The 
direct consequence of this fact is that, in our model, the effects we measure cannot be ascribed to the variation 
a(t) alone. In particular, in a continuosly varying-c theory, it is incorrect to identify the shift in the frequencies 
of a photon with the evolution of the scale factor a(t). 

As anticipated, a way to tackle this problem is to account for such additional dynamical effects by providing 
the relation existing between the evolution of the scale factor and that of the apparent scale factor. To accomplish 
this task, the first step is to calculate d^/ g rr / dx° in terms of a(t). What we obtain is 

d \Zg^ _ d = d a = a , , 

da; c(t) dt + 1 dc ° U (c + tc) dt ° (c + tc) ' ' 



This result will provide us with a connection between the observed and the actual dynamics of our universe. 
Let us denote this quantity with a — a(t) just to have a concise notation^] for the last term in equation Q. We 
may then rewrite the last equation as 

a = 7 rr (8) 

In what follows we will show that this variable works like a sort of bridge between a(t), the dynamics of the 
universe, and A(t), which is what is actually obtained when estimating a(t) from measured frequency shifts of 
the radiation coming from distant objects. Before we can give the exact relations which link these two variables 
and their derivatives, we need to calculate the basic tensor quantities for the present model. This is done in the 
next section. 



4 Basic tensor quantities 



The new metric in spherical coordinates for a FLRW universe with a time-dependent speed of light is: 



f[c(t) +tc(t)]" 

a 2 (t) 



1 — kr 2 









-r 2 a 2 (t) 









r 2 a 2 (t) sin 2 (6) J 



(9) 



5 We choose not to use the prime ' to avoid confusion with the established notation which denotes the derivative with respect to 
conformal time. 
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which yields the following connection coefficients with a contravariant t index: 

\ 



/d t [c + tc] 
c + tc 






at 



V 



(1 - kr 2 ) (c + tc) 






•\2 





r 2 aa 
(c + tc) 2 






r 2 aa sin 2 (6») 
(c + tc) 2 J 



(10) 



The connection coefficients with a spatial contravariant index are the same as those for the FLRW metric. In 
fact gtt depends on t only and no further changes have been made in the usual FLRW metric, thus the "spatial" 
T's are left untouched. 



The Ricci tensor is diagonal and its components are: 



Rt. 



-3 (c + tc) d 
a dt 



c + tc 



(11) 



Rr 



(1 - kr 2 ) 



2k 



a (c + tc) dt \c + tc 



aa 



(12) 



Rgg = r 2 (1 - kr 2 ) R r 



= r 2 sm 2 (9) (l-kr 2 ) R rr = sin 2 (6>) R ee 



(13) 
(14) 



The components of the Einstein tensor are: 

3 



G 



*~ a 2 



(c + tc) 2 



(15) 



Gr s~i9 s~i<p 
, r — \J n — Lr , — 



k + 



+ 



(c + tc) dt\c + tcj (c + tc) dt \c + tc 



1 d 

act dt 

1 d 



c + tc 



a dt \(c + tc) 2 



c + tc 



2a 



a 



(c + tc) dt \c + tcj (c + tc) 2 



a(c + tc) 2 dt 
i 



[log(ad)] 



a 



1 d_ 

a dt 



aa 
c + tc 



(16) 



where we have given several different expressions for G r r = G e e = G"^ as each of them turns out to be 
useful depending on the situation. In some of these expressions non- vanishing denominator(s)/argument(s) are 
assumed so their use is subject to a preliminary check of such condition(s). 
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Finally, the curvature scalar is conveniently expressed as: 



R 



6 



1 



(c + tc) dt V (c + ic) 



(17) 



5 Apparent acceleration in a flat universe 



Here we establish the existence of apparent effects that, in our model, are measured by constant-c observers. 
We discuss the relation between the actual and the apparent evolution of these scale factors. We will show that 
such a universe appears to be in an accelerated expansion (A, A > 0), when its actual motion is that of an 
accelerated contraction, i.e., when the scale factor a(t) is a decreasing function of 

As already remarked, any measurement made by observing some kind of radiation coming from another event 
in space-time is a test of the evolution of both the spatial and the time components of the metric. Therefore, 
what we effectively estimate is the apparent motion of the universe as desumed from frequency measurements 
performed by our instruments. Of course all of the existing estimations have been performed without taking 
into account a possible variation of c(t) with cosmological time. It follows that any direct comparison between 
observations and the theory are intrinsically flawed if c(t) is time-dependent. In fact, a direct comparison 
inevitably leads to a wrong guess for the behaviour of a(t). More precisely, it leads to a confusion between a(t) 
and Alt). 

What we set out to prove now is that a homogeneous and isotropic universe with a varying speed of light 
and filled with "normal" content only (that is, without exotic matter or a cosmological constant) may appear 
to be in an accelerated expanding motion, this apparent effect being caused by a variation of the speed of light. 

We provide two ways to show what we have just claimed. The first way is rather simple and makes use 
of the relation existing between the wavelength and the frequency of a given wave which, in the case of an 
electromagnetic wave propagating in vacuum is: \v = c. It consists in calculating the apparent and the real 
evolution of the wavelength starting from a given measurement of the frequency of the same wave and considering 
the variation of c with time. It is clear that this method is well founded only if we are making a comparison 
between what is recorded when observing radiation coming from a distant object and what is known about 
atomic physics and atomic spectra. The second method is more general and applies to any possible case. The 
two methods give the same results. 



5.1 Motivations for space-time flatness 

Before looking into the first proof we need an intermediate result. In fact, the starting point of both proofs is 
to require that the curvature scalar vanishes. In this paper we limit ourselves to the study of a spatially flat 
universe (k = 0)0 so that the requirement becomes 

R(k = 0) = 0. (18) 

This assumption may appear controversial so we will dedicate a few lines to discuss it in some detail. 

First we remark that, while the vanishing of the Ricci scalar is assumed, we can easily argument that the 
Ricci scalar, whatever its value, should remain constant. In fact the gravitational interaction, in its geometri- 
cal representation, is described by the dynamics of a four-dimensional manifold which, in a sense, "mimic" a 
gravitational field. The curvature of this manifold is expressed by the Riemann tensor. By contracting this 
tensor two times, one obtains the Ricci scalar which, therefore, is the scalar quantity that best represents the 
curvature of the manifold [35] . 

6 As we will soon show, this is not in contradiction with the Hubble law. 
7 The spatially curved case will be discussed in a paper under preparation. 
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In a cosmological setting, for a homogeneous and isotropic universe, the form of the metric is that of 
Robertson- Walker, i.e., that of a maximally symmetric space so that, in the FLRW as well as in our model 
(which can be seen as a VSL extension of the FLRW model) the Ricci scalar R can, in principle, depend on 
time t only, and not on the spatial variables: R = R(t). 

Next consider that General Relativity is based on the Einstenian conception that gravitation in nothing 
but the curvature of space-time, and that this curvature can be generated by mass-energy-momentum and only 
by it. Given that, by its very definition, the universe cannot interact and exchange mass-energy-momentum 
with anything else, it immediately follows that the average space-time curvature of the universe should not be 
allowed to change with time. As the space-time curvature can be expressed, for instance, by the Ricci scalar - 
which because of the cosmological principle cannot depend on the spatial variables - we arrive at the conclusion 
that the Ricci scalar must remain constant in time. 

We stress that all this discussion applies to a homogeneous and isotropic cosmological setting only. In fact, 
in different situations, the possibility of exchanging energy-momentum with an "outside" forbids to follow the 
above line of reasoning and therefore the Ricci scalar cannot be constrained if not in a cosmological ambit. 

The fact that R must remain constant does not however mean that it must vanish. The value zero has 
been chosen essentially for two reasons. First, somewhat philosophically, we may say that it seems to be a 
natural choice: an initial condition must be chosen and deciding that the universe "has come into existence" 
with no prior space-time curvature appears to be, in a sense, the most natural one. Second, more formally, if 
we want a classical universe to emerge from its quantum epoch (and from a possible inflationary stage) in a 
radiation-dominated state, then we must fix R = R RAD as a sort of matching condition. As during this stage 
the trace of the energy-momentum tensor vanishes T = = it follows that R — for any model based on 
Einstein field equations. When this is combined with the previous remarks, we are led to fix R = once and 
for all. 

We will soon see that the fact that some forms of matter, e.g. dust, do not satisfy T = = (which is 
implied by R = 0) does not represent a problem for our model. Indeed the variation of c will be seen to induce 
additional terms in the components of the energy-momentum tensor. These additional terms will in turn induce 
a solution to T = = which is different from e = 3p or from the trivial one, e = (here e is the energy 
density) , so that fixing R = does not constrain the equation of state in the model we are proposing. 

We digress for a moment to note that the assumption R = 0, when made in the standard FLRW cosmology, 
enforces a radiation-like equation of state (or that of an empty universe). In fact, -R FL r W = implies a(t) oc 
t 1 / 2 and this shows that a flat FLRW space-time seems to be compatible only with a permanently radiation- 
dominated universe. The "creation of a net amount of curvature" taking place locally at every point of the 
spacelike hypersurface of simultaneity in the standard FLRW cosmology at the transition from the radiation 
to the matter dominated epoch (during which the Ricci scalar decreases as 1/t 2 ) is hard to explain when one 
adopts the viewpoint outlined above. Since gravity couples to any kind of energy-momentum, and since the total 
amount of energy-momentum must obviously be conserved during the transition from the radiation-dominated 
ot the matter dominated-phase, it follows from the contracted field equations that the mean space-time curvature 
represented by the Ricci scalar must be conserved as well. As already said, this problem does not arise in this 
VSL model, where one can impose the vanishing of the Ricci scalar and still retain the possibility to have any 
kind of equation of state that satisfies the positive energy condition. 



5.2 Consequences of space-time flatness 

The assumption of space-time and spatial flatness, when combined, turn out to be very useful. From equation 
(jTTI) , keeping only the last term which does not contain k, we get 



R 



-6 



a 2 (c + tc) dt 



(c + tc) 







(19) 
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which provide us with a fundamental relation between the scale factor and the speed of light : 



aa 



(c + tc) 

The general solution for the scale factor to equation l|20p is 

a(t) = a 



CONSTANT. 



tc 

*0 Co 



(20) 



(21) 



and it is immediate to realize that the only elementary function fulfilling such requirement is a power law. This 
could have also been seen by noting that, for the left hand side to be constant, that is, independent from t, the 
first requirement is that both terms in the denominator must have the same dependence on t. This can be so 
if and only if c(t) evolves according to a power law. This forces the numerator to be also a power law function. 
So, both the scale factor and the speed of light must evolve according to a power law if the evolution law for 
the universal expansion has to be simple. 



Imposing for a(t) and c(t) the following forms: 



a(t) 
c(t) 



with N, M <E R, we get 



which immediately yields 



t-M 



Nt 



aa_, N 

°0 .M 
j.M 1 
l 



2N-1 



cq J (1+M)t 



M 



CONSTANT 



M = 2N - 1 



CONSTANT = 



2t C 



so that equation l(20|) becomes 



(c + tc) 2t c 



(22) 
(23) 

(24) 

(25) 
(26) 

(27) 



We remark that using ^fg^t = c(t) instead of (c + tc) does not allow to obtain the previous relations so that 
the whole of the results that follows cannot be extended to other continuosly-varying speed of light theories. It 
is also not possible to draw the conclusion that c(t) and a(t), must be power law functions. 



5.3 An intuitive argument 

Keeping the last result at hand, let us now suppose that we have measured the properties of the radiation 
coming from a set of distant objects. In particular let us assume that we have measured the frequency of a 
given emission or absorption line, compared it with what we already know from atomic physics about that 
process, and have finally deduced a dependence of that frequency on cosmological time which can be written as 
v oc t~ s with SeK. As a deduction, if thinking in terms of a constant speed of light, we get that Ait) oc A oc t s 
(let us remember that A(t) is the apparent scale factor, apparent in the sense that it is what is observed by a 
constant-c observer). 

Now let us think in terms of a varying speed of light. We have just seen that, if the scale factor behaves 
as t N , then the speed of light is proportional to i 2JV_1 , see (f23|) and (J2SJ) - From Xu — c, we then obtain that 
frequencies are proportional to t N ~ x . 
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As our measurements gave v oc t~ s , we get that N = 1 — S and not S is the correct exponent in the power 
law evolution of the wavelength and thus for the scale factor. In other words the scale factor evolves like 
but appears to be evolving like t s : 

IF a(t) OC i 1-5 THEN A(t) oc t s . (28) 

As 

A{t) > < — ► S < V S > 1, (29) 

we deduce that a positive acceleration is detected when N > 1 or N < 0. In particular the second possibility 
tells us that a contracting universe with a varying speed of light manifests itself as a universe that is undergoing 
an accelerated expansion, and this is exactly what we set out to show. 

The other solution, yielding N>1, is to be discarded for several reasons. For instance, it implies S<0, that 
is, the observation of a contracting universe which is in contradiction with the experience. Also, it implies an 
increasing speed of light, which is really unwanted as the speed of light must decrease for the horizon problem 
to be elegantly solved by a VSL mechanism. But most importantly we will not consider the possibility TV > 1 
since, as we shall see later, it will be formally excluded upon solving the field equations. Indeed, it will be shown 
that each solution to the field equations satisfies N e [—1/2, 1/2] which clearly excludes N>1. 

We remark that the result given above concerning the apparent nature of the acceleration has been achieved 
with no assumption about the content of the universe. No exotic matter and/or cosmological constant have been 
introduced. The effect is simply induced by the variation of the speed of light. 



5.4 A general proof 

The argument exposed in the preceding section is not rigorous and may not apply in some cases. There is 
however a more general way to show how and when an observed acceleratation turns out to be an apparent 
effect. We make use of (0 to calculate the evolution for a in terms of the cosmological time: 

■ h a ° t»-H- N (30) 



(c + tc) 2c Q 

The evolution of a is a way to represent the dynamics of the manifold and, in a sense, it provides a bridge 
between the real and the apparent evolution of the scale factor expressed in terms of cosmological time. In fact, 
an observer who does not take into account the variation of c(t) would perceive an evolution for the apparent 
scale factor given by 



A = ch = w t ■ (31) 



a_o_ j.N- 

The crucial point to note is that A is always positive, irrespective of the value of N. This is a very important 
property of this model: no matter whether the universe is expanding or contracting, the observed (apparent) 
motion will always be found to be that of an expansion. 

This means that we may well be living in a varying-c universe whose spacelike hypersurfaces are undergoing 
a contracting motion, but such a universe would anyhow show just the same Hubble law a standard FLRW 
universe shows. Therefore, and we stress this very important point, the fact that N < is not in contradiction 
with the experimental evidence that we observe redshifts and not blueshifts. 

The result for the apparent second derivative of the scale factor is even more interesting. We make use of 
the previous result for the apparent first derivative of the scale factor to derive the apparent acceleration: 



As this seems to be a source of confusion, we remark that A is not the time derivative of A just as A is not 
the time derivative of A. In words, A is the apparent time-derivative of the scale factor, not the time-derivative 
of the apparent scale factor: A = (a) APP ^ -^(a APP ). 
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The previous equation shows that an acceleration for the expansion of the universe is measured if and only 
if N<0 V N>1 : 

POSITIVE ACCELERATION OBSERVED N < V N>1. (33) 

Note that this is exactly the same result we have obtained with the first method. We made no use of the 
properties of a wave propagating through space-time. It is a simple computation of how modification of distances 
induced by the cosmological expansion are seen by an observer who does not take into account the variation of 
the speed of light. 

The crucial point to note is that A can be positive even when N is negative. It is precisely this result that 
allows N to be negative, thus causing a positive acceleration, without the need for A to be negative as well, 
which of course would mean to observe a contraction, in contradiction with the experimental evidence that 
we observe redshifts and not blueshifts. This is evident by looking at the summary for the evolution of the 
(apparent) derivatives of the three variables a(t), a(x°) and A(t) which is given in Table [IJ 





(Apparent) 1 st derivative 


(Apparent) 2 nd derivative 


a(t) 


Nt N -' 


> for N>0 


N(N- l)t"-' 1 


> for A<0 ViV>l 


a(z u ) 


1 + -N 

2 L 


> ViV 


4 1 


< ViV 


A{t) 


1 jjj-l 

2 1 


> ViV 


N-l ±N-1 
4N L 


> for N<0 V N>1 



Table 1: Time evolution for the three "scale factors" in terms of time. Note that a < if A< but A > VAT. 
We stress that A is not the time derivative of A just as A is not the time derivative of A. In words, A is the 
apparent time-derivative of the scale factor, not the time-derivative of the apparent scale factor. 

Hidden behind the above line of reasoning is the assumption that there is no variation of the emission 
properties with cosmological time. Indeed, this feature derives from the choice of the time variable and is 
therefore automatically guaranteed in our modelH If that was not the case, we could not make any comparison 
between atomic spectra obtained on earth laboratories, and emission/absorption lines of distant objects. This 
issue has already been discussed in section [2] but a further discussion is perhaps at order, given that the 
characteristics of the emitted radiation depend on the fine structure constant a EM and that, as mentioned in 
the introduction, some recent papers (the last one being [15]) have reported a very tiny relative variation of this 
constant. We postpone this discussion to section [7] where the consistency of the model with observations and 
its predictions are addressed. 



6 Field equations and their solutions 

We have shown that a universe that appears to be undergoing an accelerated expansion actually is, according 
to the present model, in a state of contracting motion. Yet we still have to prove that such a behaviour is a 
solution to the field equations of the model and discuss what the content of the universe has to be in order to 
obtain it. Our next task will then be that of solving the field equations. However, one runs into difficulties right 
from the start. In fact in this model an additional problem arises, that of expressing the energy-momentum 
tensor in terms of pressure and energy density. This issue must be tackled before we can even write the field 
equations. 

6.1 The energy-momentum tensor 

The problem is easily explained as follows: in standard cosmology one has T°° = po which after lowering an 
index using g 00 = c 2 , becomes T° = e = poc 2 . Here po is the mass density in the rest frame of the fluid and e is 

8 Note that the same is instead assumed in standard cosmology. 
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the related energy density. This shows that one can define the energy-momentum tensor in terms of the "usual" 
quantities (density of energy, of mass, momenta, pressure, etc.) using either contravariant or mixed indices. 

This remains true in VSL theories, even if p and e do not have the same time dependence anymore. In this 
model this arbitrariness is completely lost. To get through this ambiguity, let us use, as a guidance, the way 
the metric tensor is defined. The metric tensor is defined in terms of its covariant components because 
multiplication by differentials dx M <g> dx 11 (product of basis 1-forms) must give the line element which is a scalar. 
The other versions, if needed, are to be derived via the usual indices raising procedure. 

By the same token, the components of the energy-momentum tensor are to be defined in their contravariant 
form. This is easy to see in the simple collisionless case where the energy-momentum tensor is defined starting 
from four- velocities in the comoving frame: T^ v = po(3^(3 u . Given that four- velocities are four- vectors (defined 
with contravariant indices) the energy-momentum tensor is opportunely defined with contravariant indices as 
well. It is only by lowering an index that we can express T^ v in terms of the energy density of the fluid. In our 
model this gives: 

T° = p(c + tc) 2 = -l(c + tc) 2 = 4N 2 e. (34) 
cr 

The reason why we have to solve the field equations in their mixed-indices form is that the energy density e 
must appear in the equations in order to use the equation of state w = p/e. 

Straightforwardly extending what we have just seen to the case where a pressure is present and the energy- 
momentum tensor is defined by 

T^ = (p+^)^/3"-pg^, (35) 
we see that the energy-momentum tensor in the mixed-indices representation becomes 



/AN 2 e 








°\ 




(4N 2 e 








^ 





-p 













—we 














-p 













—we 





V o 








-p) 












—we J 



(36) 



We stress the great importance of having both N and w appearing in the mixed representation of the 
energy-momentum tensor. Without this feature, this model would have had exactly the same problem that 
affects standard cosmology and that has been already discussed: a change in the equation of state would imply 
a change in the average curvature of space-time not deriving from a change in the total energy-momentum 
content of the universe, in a flagrant contrast with the very foundations of General Relativity. 



6.2 Field equations and Bianchi identities 



If one writes the Einstein-Hilbert action for General Relativity in terms of a; and variates it, one finds that the 
Einstein field equations are not affected by the introduction of a varying speed of light. As a consequence, no 
additional terms appear in the field equations which therefore have the usual form. Using the energy-momentum 
tensor given by equation lj36|) , the quantities given in Section H] and introducing the shorthand notation 



X = X(t) = {c+ tc) 
we find that the general field equations for this model are: 



X 



IttG 



4N z e 



(37) 



(38) 
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r • 2 "i 

aa 
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a 


dl 


X 2 



8ttG 



8ttG 



■ w e 



(39) 



where e is the energy density, p is the pressure and w is the ratio p/e representing the equation of state as usual. 
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Imposing the once-contracted Bianchi identities = 0, we see that the usual conservation law for the 

energy-momentum tensor is no more satisfied: 



T M + 0. 



(40) 
(41) 



The Bianchi identities give instead: 

(k E T>i). ifi = 0, 

the "spatial" equations being trivial as usual, while the v — equation for a constant N ^ is calculated and 
conveniently rearranged into: 



w „ a , c 



G 



4iV fc -U + 3(1 + i^^" 4 c + G ^° 



(42) 



Apart from the trivial solution e(t) = 0, the previous equation has the general solution: 



Go c\t) 



4 G(t) \a(t), 

where the index zero refers to an arbitrary time as, for instance, the present timed 
If N = 0, equation l[4*Tj) becomes 



3iu feg; e — = 
a 



which is an identity given that a(t) oc t N . 
Using equation (|4*3l) the field equations become 



k + — 

X . 



= 4iV k Eo £o ( — , 



(43) 



(44) 



(45) 
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aa 




^2 ■ 
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. X 2 _ 





/ao\ 3 ( 1+ if2 
wk Eo e — 



v a / 



(46) 



where fcg is the Einstein constant computed at the initial conditions imposed upon solving equation l(42j) . its 
value being: 

k Ea = — — • (47) 



We remark that the evolution of c(t) and G(t) does not affect the right hand sides of the field equations. 
Actually G(t) is completely disappeared from the field equations (only Go, inside k Eo , is still present) so that a 
possible evolution of G(t) with time has no effect on the dynamical. 

We wish now develop further the analysis of the space-time flat solutions. We will show that indeed this 
model admits contracting solutions that can match observations without requiring the presence of any kind of 
exotic matter/energy, i.e., assuming that only radiation or pressureless matter are present. 



9 As is commonly done in standard cosmology, we have assumed that w does not vary appreciably over times where the solution 
given is used, i.e., we are studying the dust or the radiation-dominated eras. 

10 We are assuming that we can introduce a variable G(t) directly in the field equation. If the variable G(t) is introduced in the 
action then the field equations will be different form those given above. The results obtained in the paper remain valid in any case 
if one reverts to a constant G. 
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6.3 Solutions from space-time flatness 



Under the requirement of space-time flatness, the most direct way to obtain information about the solutions 
to the field equations is to take advantage of the fact that the trace of the energy-momentum tensor vanishes: 
T = G = —R = 0. This requirement immediately yields AN 2 — 3w, which, solved for N, gives 

N = ±^. (48) 

Needless to stress how important this relation is, and what are the implications of the "±". It directly gives 
the exponent for the power-law behaviour of the scale factor in terms of the state parameter w. 

Making use of this relation one obtains the following equations of motion for the scale factor and evolution 
for the speed of light: 

a(t) = a ^-J (49) 

To) • (50) 

These are to be compared to the solution one obtains for the standard FLRW universe: 

/ t \ 2/3(1+*;) 

o(t) = a f-J (51) 
c(t) = c . (52) 



Let us now analyse the properties of the solutions we have just found. The first feature one notes is that 
we recover the same solution for the expanding radiation-dominated era, a(t) oc t 1 / 2 , but also obtain the 
solution a(t) oc t^ 1 / 2 , which is the contracting solution we were looking for. As for the matter era, we find 
a(t) = constant. We remark however that the constant solution is exactly constant if and only if nothing 
but pressureless dust contributes as source into the energy-momentum tensor. Even the smallest contribution 
from any form of energy (e.g. radiation) different from the pure dust is sufficient to give a varying scale factor. 
Therefore, this static solution is just an idealized and asymptotic one. 

The crucial point to note is that the general solution given in equation (|4"8)l present us with two solutions, one 
being the reciprocal of the other. The most evident implication is that a universe containing any given mixture 
of dust and radiation (with one dominating over the other) can be either in a contracting or in a expanding 
motion, and which of this two motions is actually taking place cannot be told by estimating the content of the 
universe, because the same w is compatible with both solutions. So, for any given equation of state, one always 
has two possible solutions for the state of motion of the universe, a contracting one and an expanding one. 

On observational grounds, one could think that the motion effectively experienced by our universe is that of 
an expansion, but one would then be forgetting TableHJ which tells us that both an expansion and a contraction 
are perceived as an expansion by a constant-c observer. Indeed, the apparent first derivative of the scale factor, 
A, is positive for any N. But if both "thermodynamical" and "kinematical" observations cannot tell an expanding 
from a contracting universe, does it make any sense at all to distinguish between the two solutions? 

The answer is that there actually is a kinematical way to tell one solution from the other, but it is not 
related to the first derivative. Rather, the imprints of a contraction are to be read from the apparent second 
derivative. Looking again at Table [TJ we see that if N < or N > 1 then the apparent second derivative turns 
out to be positive, while for < N < 1 it is negative. Observations tell us that A > 0, so that our universe 
is now a mixture of dust and radiation-like matter and is undergoing an accelerated contraction powered by 
energy-pressure. That is how the puzzle of the acceleration of the universe becomes the only instrument we have 
to tell in which phase, expanding or contracting, the universe presently happens to be. 
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6.4 An eternally bouncing universe? 

It is interesting to speculate on what the global history of such a universe could have been and what its future 
may be. A possible scenario is that of a universe undergoing an eternal cycle of successive expansions and 
contractions. This is similar to what happens in the standard FLRW universe with a positive spatial curvature, 
but with the two main differences that here the universe is both space-time and spatially flat, and that the 
evolution law for the scale factor is not a cycloid but rather a power law with a varying exponent. 

To give a brief qualitative description of such possibility, let us start from a hot and dense universe where 
radiation dominates. Once the radiation era starts, the universe expands according to the well known power 
law a(t) oc i 1 / 2 , with the exponent decreasing as the expansion cools the universe until pressureless matter 
dominates over radiation. In this phase the expansion becomes very slow (the exponent tends to zero). 

Since pressure is ever decreasing, general relativity approaches a Newtonian behaviour, where mass density 
generates what looks like as an attractive force. For some reason (e.g. because of substantial deviations from 
a pure homogeneous and isotropic scenario), the expansion might eventually stop and an almost-classical con- 
tracting motion could begin. Such motion is described by the same power law but with an exponent of opposite 
sign. According to this almost-Newtonian behaviour, the contraction accelerates more and more (exponent in 
the power law increasing more and more as the universe contracts), until the pressure regains his role in the 
game. 

Eventually, the universe enters a new radiation-dominated era and, according to our solution, the exponent 
N stabilizes to a value equal to —1/2. The now constant but steady contraction will finally cause the matter 
density and the pressure to reach such high values that quantum effects become important and classical general 
relativity loses its fundamental character. Possibly, during this quantum era, the universe may experience a 
sort of Big Crunch followed by a new Big Bang which makes the cycle to start again. 

6.5 Evolution of a(t) and c(t) and some related speculations 

Given that g tt = X 2 = 4iV 2 c 2 (i), that equation (08]) implies 

tcG[0;l/3] — ► Ne [-1/2; 1/2] (53) 

and that 

c(t) = c 

we find that both g tt and c(t) are non- increasing functions of t for any TV e [— 1/2; 1/2]. In Table [2] a summary 
of the asymptotic evolutions of c{t) and a(t) in both the contracting and the expanding phases is given. 





a(t) 


c(t) 




Expansion 


Contraction 


Expansion 


Contraction 


w = 










w = 1/3 


f l/2 


f -l/2 


t° 


t- 2 



Table 2: Evolution of a(t) and c(t) for w = 0, 1/3. 

An interesting feature is that, for the radiation-dominated era in the expansion phase, we recover the 
standard FLRW cosmology: a(t) evolves as t 1 / 2 , 47V 2 = 1 and c(i) is constant. The FLRW model matches 
then with our model during a small but important part of the history of the universe. Most importantly, this 
guarantees that the description of the post-inflationary early phases of the universe in our model is precisely 
the same of the classical FLRW description. For instance, results concerning the BBN are not affected by a 
varying-c as c(t) does not vary at those high energies in the present model. Note that this is the opposite of 




17 



what happens in other popular VSL models, where the speed of light often varies only at very high energies 
(during a phase transition). 

Apart from this brief period where c(t) is asymptotically constant, the speed of light is a decreasing function 
of time. c(t) varies according to a power law with the exponent ranging from the idealized value of 0, during the 
radiation-dominated expansion phase, to a maximum value of -2, during the radiation-dominated contraction 
phase. This guarantees that any given part of the universe was causally connected to any other different part 
at some moment in the past. No matter how far two points in space are nowadays, one can always find an early 
enough time when these two points were causally connected. In this model there is no horizon problem. 

Another curiosity we wish to address is that for N =±1/2 - that is, during the radiation-dominated eras - 

we have x(i) = c(t) so that one recovers the "usual" line element ds 2 = c(t) 2 dt 2 — a(t) 2 dr (usual in the sense 
that c(t) appears instead of x(i)). This suggests the possible existence of an interconnection between the fact 
that photons (massless particles) travel along null geodesies, the fact that they travel precisely at a speed equal 
to c(t) and the fact that the pressure of radiation is equal to 1/3 of its energy density. In fact, if any of the 
previous things is changed, then the other two will change as well: all three are concomitantly true or false. 
Such result suggests that some properties of massless particles may have a cosmological origin. 

The dust-dominated era is not less intriguing: given that x = 4iV 2 c(i) 2 — > when N — > 0, we have a sort of 
"cosmological indication" that massive particles have smaller speeds with respect to the locally measured speed 
of light. It is tempting to use naive kinetic theory to transform the line element. From the usual formula: 

P=^("rms) 2 (54) 

we obtain 

w = ^(%ms) 2 (55) 
that, once substituted into g u — 4iV 2 c 2 (i), making use of iV 2 = 3io/4, gives 

9tt = (vrms) 2 (56) 

so that the line element becomes 

d S 2 ^ 2 MS (i)dt 2 -a(t) 2 dr (57) 

which seems to suggest that particles of any mass, depending on the idealized content we attribute to the 
universe, can define worldlines of null length. 

Of course we are abusing here of the concept of v RUS by assuming that any particle in the swarm has the 
same velocity. If we however were to take all the particles to have the same velocity (a very highly idealized 
application of the cosmological principle), we could drop the RMS index and what has been conjectured above 
would become a little bit more precise. 

This discussion is of course highly speculative and would only make sense in a hypothetical universe where 
one and only one sort of fluid were present and where all the particles of that kind were to travel at the 
same speed. Consequently, the remarks made above are in no way not related to our universe. However, some 
interesting effects could take place even when a more realistic case is considered and this is why I have mentioned 
these reflections here. 



6.6 Initial conditions 



We now solve the field equations directly to look for constraints on initial conditions. We recast equation ([27|l 
in the form 



a Qq 1 



X 2£ c a 



(58) 
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and use it in equation l|45p with k 



to obtain 

2 







2i c 



(59) 



We see that a(t) drops out form this equation if and only if 3w — AN 2 , which is precisely the solution to 
equation l[4*8|) . In this case simplifications yield 

3 = (4Nt c ) 2 k E0 e 



that can be recasted, for instance, into 



1 



4JVc V k Ea e 



(60) 



Using the accepted numerical values referred to the present time for the quantities involved (assuming that 
such estimations are not excessively model dependent, which is far from being granted) one obtains 



1q26^27 

t ~ — — — sec 

the interpretation of which cannot of course be straightforward for a variety of reasons. 



(61) 



For instance, with gtt depending on time, the cosmological clock is seen to tick slower or faster depending 
on the epoch the observer happens to live in 11 !. This is reflected by the presence of N in the denominator of 
the previous equation. In fact, for N — > (matter era), to — > oo in agreement with the fact that in this model 
the cosmological clock is seen to tick slower during the matter era (g u — > when N — ► Q). 

Moreover, g tt contains the speed of light c(t), which has been showed to be a non- increasing function of time. 
This implies that, apart form the effect discussed above, the cosmological clock "ticked faster in the past" so 
that, in a sense, to is expected to be much smaller than the value given in equation (|6"Tj) From this qualitative 
discussion it emerges that the value given above cannot be directly compared to other estimations of the age of 
the universe based upon a constant speed of light (and a single-cycle universe) . 

Note that in the case 3u> ^ 47V 2 , equation (|59|) reduces to 



«o 
a 



{ANt c ) 2 k Eo e 



AN* 
w—4N 2 



(62) 



which implies that a(t) = CONSTANT and therefore a(t) = ao so that again equation <f60|) holds. However, 
equation ([27| . that we have used to simplify l[59|) . follows from the hypothesis of space-time flatness (R = 0), 
that is from 3w = AN 2 , so that one could hardly give any meaning to the 3w ^ AN 2 case. 

We turn our attention to the second field equation (|46|) . Using the second of the expressions given for G r r , 
we get 



1 




1 d 


f aa\ 2 


1 


^2 


k -\ 


aa dt 


\ X J 





/a \3(i+- 
-w k Eo £o — 



r) 



\ a / 



(63) 



where, by equation {j7?\ with k = 0, the first and the second terms on the left hand side vanish while the last 
term is constant. Under these conditions one can recast the second field equation into 



zl (jLX 

\2toCo) 



/ao\ 3 ( 1+ i#7) 
-wk Eo e , 



(64) 



11 For instance, if an observer at rest in a radiation-dominated area of the universe were able to observe the clock of another 
observer at rest in the same universe, but in a far enough portion of it which happens to be dust-dominated, then he would observe 
that clock ticking slower with respect to his own clock. 
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which both in the case 3w = AN 2 and 3w ^ AN 2 gives the same result - equation ([60)) - obtained from the first 
field equation. 

The fact that the field equations do not give additional information about the dynamics was to be expected 
as we have applied two constraints - space-time flatness and the Bianchi identities - on a system with two 
degrees of freedom: the two scale factors a(t) and I n such situations both the field equations become 

"initial conditions" equations and, as is to be expected, both give the same result. 



7 Consistency checks and predictions 

The new VSL model presented in this paper is capable of solving the puzzle of the acceleration of the universe 
but some of its aspects, such as the presence of dust and radiation only or that the acceleration is only an 
apparent effect must surely result somewhat unfamiliar and one could wonder if the model is compatible with 
the relevant observative evidences. While a complete data analysis related to the most important observative 
constraints is evidently beyond the scope of the paper, some consistency checks are perhaps at order, at least 
to be sure that no contradictions, either intrinsic or related to observations, could arise. We do not digress on 
the SN data as it is evident that the presence of an apparent acceleration makes our model compatible with 
that data set. On the contrary, there are two issues which are to be addressed with care: the critical density 
for the model, which is an internal consistency check, and, as anticipated at the end of section the relation 
with the variation of a EM , which can be considered the only observative evidence which is somewhat related to 
a varying-c. 



7.1 Critical density 

We have already emphasized that only dust and radiation are needed for our model to exhibit an (apparent) 
accelerated expansion. No dark matteJil, let alone dark energy, are needed to explain the acceleration of 
the universal expansion, which was the main issue we wished to address. However, in our analysis we have 
also assumed that the universe is spatially flat and this implies that the energy-mass density must be critical 
according to equation (|38| (with k set equal to zero) that we re-express for convenience as 

PcR = 24^ ■ (65) 

This critical density must be achieved without invoking exotic forms of energy as these have not been 
introduced into the model. At first this may seem to be a fatal shortcoming for the model as the average 
baryonic plus relativistic matter-energy density of the universe is observed to be around the 4% of the standard- 
cosmology critical density = 3H 2 /8irG implied by the standard Friedmann equation. Things worsen if 
one considers that, for a dust-dominated universe, p CR given by equation ([65| is greater than p^ L r W - Indeed one 
finds that 



From 



follows 



I v ±V3w/2 , x 2/3(1+11)) 



0(*)tHIB MODEL — a ( —J ! O^FLRW — fl \J~ j (^7) 

H 2 3u>(w + l) 2 , , 

- (68) 



H 2 16 

FLRW u 



and therefore 



KW + l) 2 FLRW ^ 

16w PcR ' [ ' 



12 Dark matter is not excluded as a further contribution to the energy-momentum tensor for the model. 
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The last relation shows that p CR = p£ R RW f° r a universe filled with radiation 13 ! but p CR > p£ R RW when dust 
dominates. After a closer inspection however, it becomes clear that things are a little more complicated. In 
fact, the critical density is calculated from the estimations made on the Hubble parameter H = a/a which, 
however, has an apparent counterpart H APP = A/A whose interpretation should be obvious by now: it is H APP 
and not H that is inferred from measurements and is therefore incorrectly used to determine what indeed is an 
apparent critical density. By means of H, we must now calculate the real critical density which is what is to 
be compared with other (redshift-independent) estimations of the matter density to check whether this model 
is consistent with observations or if it necessitates a generalization to include other forms of matter-energy in 
order to cope with observative data. 

A convenient way to follow is to calculate the ratio between the critical density as it appears from observations 
and the real critical density, that is, the matter density which is actually needed to have a spatially flat universe. 
The computation is straightforward and yields, for the present timd* 4 ]. 

PcR 



APP 

r cr 



3w (70) 



whose order of magnitude may range from 10 _1 to 10~ 4 depending on several issues as, for instance, neutrino 
contribution to the present-day energy density and pressure. In other words, the true critical density is a 
factor 10 _1 to 10~ 4 smaller than the estimated value. This broad range is in a qualitative agreement with the 
above-mentioned 4% discrepancy. 

Summing up, we have found that a determination of the critical density made by estimating the Hubble 
parameter leads to a large overestimation during the matter era. Much less mass-energy is needed to have a 
spatially flat universe with respect to that suggested by the measure of the Hubble parameter. This finding 
makes the model internally consistent and provides an elegant solution to the problem of the missing mass / energy 
with no need to resort to hypothetical forms of matter/energy. We have constructed a spatially-flat, "normally"- 
populated, apparently-accelerating universe. 



7.2 Present-time density contribution from relativistic matter 

Equation <[70|) predicts the existence of a contribution from relativistic matter to the present-time matter density 
and allows to obtain a rough prediction of its value. This contribution is necessary for the model to be self- 
consistent. We assume that there is no need to introduce any dark matter into the model and denote with 
f2 N the observed fraction of the "normal" matter density p with respect to the actual critical density p CR - The 
currently accepted value for fi N is f2 N = 0.04. If the universe is spatially flat the observed matter density must 
be equal to the apparent critical density, so that 

pobs = pr (71) 

must hold. From equation {TO} it follows then that it must be 

3w = ft N (72) 

whence 

Wtoday ~ 0.013. (73) 

Given that for relativistic matter p REL = e REL /3, we find 

R = J 
e 3e T 



P e REL /_.n 

w = - — (74) 



which, combined with equation ([72} , gives 



N (75) 



13 This further confirms that, as we had found previously, the FLRW model and this model are indistinguishable during the 
expanding radiation-dominated era. 



14 For arbitrary times the solution is 3ui(to/t) 2 ' v ^™ 
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The last equation tells us that the present-time contribution of relativistic matter to the total matter density 
is predicted by the model to be of order 4%- Such contribution must originate from "normal" (i.e. well known 
and observed) relativistic matter and the optimal candidate for this role is the neutrino. Given the uncertainty 
about neutrino masses this 4% contribution is well within the observative bounds. Actually, the previous result 
may be of help in determining neutrino masses and number densities from cosmological measurements but a 
more extensive analysis is beyond the scope of this paper. 



7.3 Connection with the a EM results 



As anticipated at the end of section the existence of the apparent effects discussed throughout the paper relies 
on the assumption that the a EM results do not have a large impact on the evolution of emission frequencies. 
In this section we discuss this issue and show that our results are correct whatever the evolution of a EM . We 
also discuss whether the variation of a EM can be explained in terms of the variation of c only. As we are about 
to see, this problem is in general ill-posed but if we specialize to the units/frame we have adopted, we find 
that the variation of a EM induced by that of c is much larger that the observed one. This finding forces us to 
conclude that in our case e and/or h must also vary. Note that this fact has no impact on the results that we 
have obtained so far as these are not affected by the evolution of h and/or e. 

We first point out that, contrary to what one may be led to think at a first glance, the variations of c(t) 
and of a EM do not necessarily imply each other. Indeed one may think that, since a EM = e 2 /he, if c varies a 
variation of a EM is to be expected as well. This first impression is wrong. As recalled in the introduction, a EM is 
a dimensionless quantity so that its variation (or non-variation) is independent from the units chosen. On the 
contrary, whatever the behaviour of a EM , one may choose the units/frame so that c (and/or e and/or h) vary 
or remain constant. 



Given that a EM oc e 2 /hc, we have 



a EM = 2 £ _ c _ h 
c h 



which, if the dependence on time (or redshift) is linear, becomes 

Aa EM = 2 Ae _ Ac _ Ah 
a EM e c h 

From this relation one sees that a EM may remain constant and nonetheless c may be time dependent, if its 
variation is compensated by that of e and/or of h. Vice versa, c can be set to a constant by an appropriate 
choice of units/frame while a EM could vary as a consequence of the variation of e and/or of h. This remains 
true in general, even if the dependence is not linear. It is because of these arguments that we affirmed that 
the question whether the variation of a vf , can be explained only in terms of the variation of c is ill-posed: it 
all boils down to the choice of units /frame 16 !. Specifying some units/frame gives a meaning to the issue and we 
will verify what happens in the units/frame adopted here at the end of this section. 

A somewhat related issue is the influence of a varying-a EM on the results that we have obtained so far: can 
a varying-a EM affect our findings in some way? As explained in section [21 we have adopted the frame in which 
the emission/absorption frequencies, as seen by an observer at rest with the cosmological fluid, are the same at 
any cosmic time. It is not difficult to see that this remains strictly true even if a EM is evolving with cosmological 
time. This is because the emission/absorption frequencies are dimensional quantities and have the dimension 
of the inverse of time. It follows that, even if a EM is a factor appearing in the theoretical expression for the 
emission/absorption frequency, it must appear multiplied by an appropriate dimensional factor in such a way 
to recover the correct dimension. Such a combination will therefore be dimensional and as such will depend 
on the units/frame adopted: even if a EM is present, one will always be able to fix the overall expression for 
the emission/absorption frequency equal to a constant. This is precisely what we have done upon choosing the 



16 Of course, if a EM is constant, a variation of c must be compensated by a variation of e and/or h. 

16 Some may consider this argument self-evident but there are still physicists who ask whether the observed evolution of a EM is 
compatible with this or that model. They do not realize that the point is meaningless. 
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frame adopted in the paper. The crucial consequence of these arguments is that any evolution or non-evolution 
of ct em is compatible with our model. Alpha results alone can in no way confirm or falsify it. 

We next examine what is the variation of a EM induced by that of c{t) obtained in the paper. We assume for 
the moment that h and e are constant and will check this assumption against the a EM results. We denote the 
said quantity with (Aa EM /a EM ) c : 



Ac 



INDUCED BY c(t). (78) 



Assuming that the variation is almost linear would be too rough an approximation so that we do no use 
equation ifTTj) and calculate the evolution directly. Observative results are given in terms of Aa EU /a EH which is 
defined as 

Aa EM = q EM (z) - q 
a EM ctQ 

where ao — 1/137 is the present-time value of the fine-structure constant while a EM (z) is the value it assumes 
at redshift z. With e and h constant a EM oc 1/c so that we obtain 

Aa EM \ c 



Using equation ([501 we find thaO 

/A \ / J- \ i+vStU 

1 (81) 



to 

In equation l|79p z is the observed z and therefore it is an apparent quantity which we will appropriately denote 
with z APP and which is obviously to be defined^ as 

1 + z AP p = ^T^y- ( 82 ) 

We then need a link between z APP and t. From equation ((2Sj) , recalling that N = 1 — S, we find that A(t) oc t 1 ^ N 
and therefore ^_ 

1 + 2 APP = (J^j ' (83) 

where we have used the contracting solution for equation l|48p . Inverting the previous relation and using equation 



81]) we conclude that 



, = (l + z APP ) - 1 (84) 

which can also be expressed as 

a E u{z A pp) = ct (1 + z APP ) . (85) 

This is the variation of a EM induced by the variation of c(t) obtained in the paper. 

Note that we have assumed throughout that w is constant over the time interval where the measuraments 
have been performed. This approximation is justified as follows. Using equation ((73|) . we find that 



_ 9 i + V3^ TODAY ; ^_ l 1 (86) 

2 + v3w TODAY 

Going backward in time in our contracting universe means approaching the asymptotic value w MAT = for 
which one gets 

_ 2 i + Vgw = _ L (87) 

2 + V3ti>MAT 



17 In equation l|50|l we are taking the minus-sign solution because in the present time the universe is undergoing a contracting 
motion. 

18 Recall that what instruments measure is a frequency shift, see discussion in section [531 
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It then follows that the assumption w=CONSTANT is well founded. Given that in this section we have used the 
contracting solution, the results presented are only valid as far as one looks back no further than the start of 
the contracting phase. 

If for z APP we use the mean of the redshifts used in [15], that is z APP ~ 1.75, and again equation l|73p . we see 
that the variation of c(t) alone would imply a much greater variation of a EM with respect to the observation: 



(2.75)" 11 - 1 ~ -0.67 (88) 



or 

a EM (z APP = 1.75) ~ 0.33 ao . (89) 
Comparing this result with the observed value [T5] 



(-0.57 ± 0.10) • 10" 5 (90) 



we conclude that, in the units/frame adopted, h and/or e must also depend on time. We stress that this fact 
has no impact on the cosmological results we have obtained in the paper as none of them is affected on the 
evolution of h and/or e. 



8 Summary 

In this paper we have developed a new varying speed of light model and have attempted to find an explanation 
for the acceleration of the cosmological expansion which does not have to resort to the introduction of exotic 
(and so far unobserved) forms of energy or of a cosmological constant. 

The model is based on the reconsideration of the most general form of the metric when the speed of light is 
allowed to vary. We have found that, if we want to ascribe the observed redshifts to and only to the dynamics 
of the space, it is convenient to work with a particular time dependent gu- Appealing to general covariance to 
perform a coordinate transformation and obtain a metric with a constant speed of light, while of course possible, 
has the drawback of destroying the above-mentioned property. On the contrary, working with the units /frame 
employed in the paper ensures that any given emission / absorption frequency is the same at any given epoch of 
the past and today, which clearly is the prerequisite to directly compare frequency shifts and to intepret them 
as purely cosmological in origin. 

We have argumented that the space-time curvature of the universe, here represented by the Ricci scalar, 
cannot have changed during its history. Naively, this can be justified because, by definition, the universe has 
nothing "external" to interact with, and therefore cannot exchange energy-momentum. As a particular case, we 
have assumed that the Ricci scalar vanishes. Our choice is dictated by the fact that the classical universe is 
believed to emerge from the quantum universe in a radiation-dominated state so that the trace of the energy- 
momentum tensor, and therefore the Ricci scalar, must vanish, whatever the model. This is equivalent to 
saying that, in a sense, the classical universe has come into existence with no net space-time curvature. It is 
to be noted that the same requirement cannot be made in the framework of standard FLRW cosmology as it 
would correspond to imposing the equation of state for radiation at any time, which would clearly contradict 
the observations. In our model, instead, this is not a problem because fixing the Ricci scalar does not fix the 
equation of state. 

We have shown the existence of the possibility that the observed acceleration of the cosmic expansion is an 
apparent effect induced by the observer not taking the varying speed of light into account. This effect arises 
because instruments measure frequencies while the dynamics of the universe is represented by the evolution of 
wavelengths and the former does not evolve as the reciprocal of the latter if the speed of light varies with time. 

Upon deriving and solving the new field equations for the spatially-flat case we have found that both a(t) 
and c(t) evolve according to a power law. For any admissible matter content - in our case for any w € [0, 1/3] 
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where w is the usual equation of state parameter w = p/e - the new field equations have both an expanding 
and a contracting solution and the two solutions are one the reciprocal of the other. More precisely, we have 
found that the scale factor evolves according to a(t) oc t N where N— ±^/3w/2. 

At a first glance, a contracting solution seems to contradict the fact that we observe redshifts and not 
blueshifts but it is not so because at fixed w the model predicts that an observer who is not taking the 
varying speed of light into account will observe the frequencies redshifted exactly by the same amount for both 
the expanding and the contracting solutions. This implies that the contracting solution is not ruled out by 
observations but, as a consequence, it also means that the fact that redshifts and not blueshifts are observed is 
not sufficient to discriminate between the two possibilities. Fortunately, the sign of the apparent acceleration is 
found to be different for the two mentioned states of motion. Thanks to this fact, we have been able to conclude 
that the observed acceleration implies that the universe is presently undergoing a contracting motion. That is 
how one of the most tantalizing puzzles of modern cosmology, the accelerating universe, might turn out to be 
the only way we have to discriminate between the two possibilities and determine the actual state of motion for 
the universe. 

Some consistency checks have been performed. We have found that the new critical density is inside the 
observative bounds for dust plus relativistic matter density, and therefore there is no need to add hypothetic 
forms of matter/energy. We have also verified that in the units/frame adopted, the evolution of c(t) implies a 
variation of a EM which is much larger than the observed one. Contrary to the common belief, we have shown 
that this does not represent a problem but just implies that, in the units/frame adopted, either e or h or both 
must evolve with cosmological time together with c(t). As e and h have no role in our computations this fact 
in no way affects our results. 

We remark that our findings have been obtained by considering the presence of dust and radiative matter 
only. No exotic forms of dark matter / energy or even a cosmological constant are to be introduced for the model 
to reproduce the observed evolution. This appears to be a very important feature as, despite the enormous 
efforts and investments, none of the mentioned forms of matter have been directly observed so far, nor has the 
cosmological constant received a satisfactorily theoretical explanation. 

Now we have a VSL model which describes a homogenous, isotropic, spatially-flat, "normally'-populated and 
apparently-accelerating universe. Such a model is clearly interesting by itself but there can be little doubt that 
it can be improved in a variety of ways. The cases of a spatially-curved universe and that of a double-component 
(dust and radiation) fluid are its natural extensions and we will investigate such possibilities in forthcoming 
research work. 
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A Proportionality between the scale factor and wavelengths 

Doubts have been posed as to whether the scale factor is proportional to wavelengths in the model proposed as 
it is in standard cosmology. Even at a conceptual level these doubts sound unfounded as the dynamics of the 
scale factor is nothing but the stretching of space. This is precisely what causes the stretching of wavelengths 
so the two quantities evolve in the same way just because, in a sense, they are the same thing. At any rate, in 
this appendix we provide a formal way to see this. 

Consider the travel of a light ray from a distant object to an Earth-based instrument and two successive 
crests denoted with 1 and 2. The space separation between the two crests is one wavelength and the difference 
of emission times is the inverse of the emission frequency. The relation linking the two quantities at the time 
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of emission as well as at the time of reception is X(t)v(t) — c(t). Instant by instant a given light ray will follow 
a null geodesic for the metric given by equation (0. If we consider the radial path from the emitting object to 
the receiver on Earth, we can fix the angular coordinates to obtain 

(c + tcfdt 2 - a 2 dr 2 = 0. (91) 

where a(t) and c(t) depend on time. Integrating equation (|9Tj) one gets: 

* :c>1 c + tc 



^rec,2 q \ -j-Q PVreg, 2 

ly-M = dr = r EM,2 - »-REC : 2 ■ (93) 

Assuming that both the emitting object and Earth are at rest with respect to the cosmological fluid, it 
follows that the two radial comoving coordinates for the emission of crest 1 and crest 2 are equal, and the same 
is true for reception. Subtracting equation |93|) from equation l|92p one gets therefore 



(*) 



^REC, 1 

dt = I dr = r BM ,i - r RECi i (92) 



dt - I — —dt = (94) 



EM. 1 



ft 



(*) «(*) 



or 

c*em.2 c _|_ rt mo ,2 c + tc 



s dt - —r^-dt = 0. (95) 

Using equations l(48|) . l(49|) . I|50p and integrating we get 

C ° (fN _ f N \ _ c (,N _ ,N \ (q(f) 

.2V-1 V REC, 2 REC, 1 / — .N — l V EM, 2 ''EM,!/ - V JU J 

ftot ao J o 

Rewriting t^ EC 2 as 

.2JV-1 

AT _ r REC,2 

<-rec,2 — T]v l REC,2 

r REC,2 

and the same for t^ EC x , t EU 2 and t EU x , and using again equations (J48j) , (|49| and lf50|) , we obtain 



ci\ / d\ / c< \ ( ct 



a I REC,2 V a 



(98) 



Given that 



we get 



Using 



Crec,2 — Crec.I — C REC ; c em,2 — Cem.I — C E m (^^) 

&REC,2 = ^REC,1 = &REC) a EM,2 = ^EM.l = ^EM (100) 

( — ) (^REC,2 — £rEC,i) = ( — ) (*EM,2 ~ £eM,i)- (101) 
V ft / REC V ft / EM 



Crec(^rec,2 ^rec,i) — Cr EC /z"rec — A RE c (102) 

and the same for emission, we finally obtain 



that is, A and a are proportional as we set out to prove. 



Arec A em (iq3) 
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